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The incompatibility of partial d occupation on the perovskite £?-site with the standard charge 
transfer mechanism for ferroelectricity has been a central paradigm in multiferroics research. Nev- 
ertheless, it was recently shown by density functional theory calculations that CaMnOa exhibits 
a polar instability that even dominates over the octahedral tilting for slightly enlarged unit cell 
volume. Here, we present similar calculations for LaCrC>3, which has the same d 3 B-site electron 
configuration as CaMnOs. We find that LaCrOa exhibits a very similar, albeit much weaker, polar 
instability as CaMnOs. In addition, while the Born effective charge (BEC) of the Mn 4+ cation in 
CaMnOs is highly anomalous, the BEC of Cr 3+ in LaCrC>3 is only slightly enhanced. By decom- 
posing the BECs into contributions of individual Wannier functions we show that the ferroelectric 
instabilities in both systems can be understood in terms of charge transfer between TM d and O p 
states, analogously to the standard d° perovskite ferroelectrics. 



I. INTRODUCTION 

The relative scarcity of multifcrroic materials has often 
been explained by a chemical incompatibility between the 
factors that promote ferroelectric ity compared to the fac- 
tors that promote magnetic order. - - Even though several 
mechanisms for ferroelectricity that are compatible with 
the simultaneous presence of magnetic cations have been 
identified recently (see e.g. Refs. [3H5]), the most com- 
mon mechanism for driving polar displacements in typi- 
cal perovskite ferroelectrics seems to require a completely 
unoccupied d shell of the transition metal (TM) cations 
on the perovskite i3-site. This statement is to some ex- 
tent based on the observation, that essentially all known 
perovskite ferroelectrics (e.g. BaTiC>3, KNbC>3, PbTi03, 
Pb(Zr,Ti)03, etc.) contain formal d° TM cations on the 
B-site. On the other hand a partial filling of the elec- 
tronic d states is required in order to create a magnetic 
moment P and the resulting incompatibility has been a 
central paradigm in multiferroics research over the past 
few years. 

The driving force behind the ferroelectric distortion 
in the d° perovskites, such as e.g. BaTi0 3 or KNb0 3 , 
is related to hybridisation between the filled ox ygen 2p 
states and the empty d states of the TM cationPEEl A 
ferroelectric displacement reduces the distance between 
the TM cation and one or more of the surrounding oxy- 
gen anions, and thus strengthens the corresponding co- 
valent bonds, while slightly weakening the bonds to the 
other surrounding oxygen anions where the correspond- 
ing bond distance is increased. This leads to an overall 
gain in covalent bond energy, which however is opposed 
by repulsive electro-static forces. The gain in hybridis- 
ation energy is maximal if the antibonding states with 
predominant TM d character are completely empty, and 
is zero if these antibonding states are completely filled. 

Nevertheless, it has been shown recently that cubic 
perovskite CaMn03 exhibits a similar ferroelectric in- 
stability which, even though it is rather weak at am- 
bient conditions, can be significantly enhanced by ap- 
plying negative pressure)^ At the equilibrium lattice con- 



stant the ferroelectric instability is suppressed by a much 
stronger antiferrodistortive instability, i.e. a collective 
tilting of the octahedral network, which leads to the ob- 
served P&nm-distorted perovskite structure of CaMnC>3. 
This antiferrodistortive mode is rather independent of 
volume, so that the polar distortion becomes dominant 
for slightly increased lattice constant. Very similar be- 
haviour has a lso be en reported subsequently for SrMn03 
and BaMnOs! 10 ! 11 ! These calculations thus predict a new 
class of multiferroics, which could be synthesised for ex- 
ample as epitaxial thin films, where strain can stabilise 
the perovskite structure with enlarged lattice constant. 
However, the question of what is the driving force behind 
the ferroelectric instability in the corresponding systems, 
in particular whether it is related to the above-described 
charge transfer mechanism, has not been addressed in 
these previous studies. 

Within a predominantly cubic crystal field, the d 3 case 
of CaMnOs can also be interpreted as e° electron con- 
figuration, analogously to the d° configuration found in 
most perovskite ferroelectrics. In this case the empty e g 
states can in principle provide the increase in bond en- 
ergy required for the ferroelectric instability, whereas the 
partial filling of the ti g states can create a local mag- 
netic moment. While the resulting gain in bond en- 
ergy is probably smaller than for the case of a d° TM 
cation, it is not obvious that this gain is necessarily al- 
ways sm aller th an the opposing effect of the short range 
repulsion! 2 * 4 * 11 ! In fact, exactly this scenario has already 
been discussed in Ref. [21 and arguments were given that 
in addition to the usual repulsive forces, the Hund's rule 
coupling would further disfavour the ferroelectric insta- 
bility in these cases. 

Here we address the question of what is the driving 
force behind the ferroelectric instability in CaMnOs and 
related systems. We clarify whether the ferroelectric dis- 
tortion in these materials is driven by charge transfer be- 
tween the TM cation and its surrounding oxygen anions, 
similar to the conventional d° perovskite ferroelectrics. 
In order to investigate whether this effect is a peculiarity 
that is specific to the Mn 4+ cation, or whether a simi- 
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lar polar instability can also be observed in other mag- 
netic perovskites, we compare the case of CaMnC>3 with 
LaCrC>3. Under equilibrium conditions LaCrC>3 exhibits 
a P6nm-distorted perovskite structure and G-type anti- 
ferromagnetic order similar to CaMn03. Furthermore, 
the Cr 3+ cation on the perovskite B site has a d 3 electron 
configuration that is iso-electronic to the Mn 4+ cation in 
CaMn0 3 . 

We calculate phonon frequencies and eigenmodes at 
the T-point and at selected zone-boundary wave-vectors 
for LaCrC>3 in the ideal perovskite structure at different 
volumes, and we then calculate and compare Born effec- 
tive charges for LaCr03 , CaMnC>3 , and the nonmagnetic 
ferroelectric BaTi03. Furthermore, we analyse the driv- 
ing force for ferroelectric displacements in these three 
systems by decomposing the Born effective charges into 
contributions of individual Wannier functions. 



II. COMPUTATIONAL DETAILS 

All results presented in this report are obtained using 
the Quantum ESPRESSO package, employing a pla ne 
wave basis set and ultrasoft pseudopotentials! 13 * 14 ! A 
plane-wave kinetic energy cutoff of 35 Ry (420 Ry) is 
used for the expansion of the wave- functions (charge den- 
sity). The 3s and 3p semi-core states of Cr, Mn, and Ca, 
as well as the 5s and 5p states of La, are included in 
the valence. Different k-point meshes are used to calcu- 
late the different quantities presented in this work, with 
convergence verified in all cases. Calculations are per- 
formed using both the Generalised Gradient Approxima- 
tion (GGA) and the GGA +t7 ap proach with U cS — 4 eV 
applied to the TM d statesPDl 



III. RESULTS 
A. Phonon calculations 

Similar to Ref. [S] we first identify potential phonon 
instabilities in cubic LaCr03 by calculating eigenfre- 
quencies of all zone-centre and various selected zone- 
boundary modes at different lattice constants. We in- 
vestigate variations of ±3% around an average lattice 
constant a = 3.89 A, which corresponds to the same 
volume per formula unit as in the experimentally ob- 
served Pbnm structure of LaCrOsP^ The relevant zone- 
boundary modes are selected by decomposing the struc- 
tural distortion between ideal Pm3m symmetry and the 
experimental Pbnm structure into symmetry adapted 
modes, and then identifying the main components in this 
decomposition. The ISODISPLACE utility is used for 
the mode decomposition.^ The dynamical matrix is ob- 
tained from the calculated forces created by small finite 
displacements of the individual ions. 

We find two strongly unstable antiferrodistortive zone- 
boundary modes (R§ and ), which are responsible for 
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FIG. 1: Calculated eigenvalues to of the dynamical matrix of 
cubic LaCrOa for selected modes as function of cubic lattice 
parameter a using U c s = 4 eV. The upper panel corresponds 
to G-type antiferromagnetic order (G-AFM), the lower panel 
to the ferromagnetic case (FM). Negative co 2 indicates a struc- 
tural instability. 



the experimentally observed Pbnm ground state struc- 
ture, and a soft polar (T7) mode, which is unstable for 
larger lattice constants.^ The corresponding results are 
presented in Fig.[T]for both the ground state G-type anti- 
ferromagnetic order and the ferromagnetic case. Similar 
to CaMnOs, the antiferrodistortive modes are rather in- 
sensitive to volume, whereas the polar mode is strongly 
volume-dependent. It can also be seen that ferromag- 
netic order leads to a further destabilisation of all modes. 
This is a result of the fact that the antiferromagnetic 
superexchangc interaction is strongest for an ideal 180° 
TM-O-TM bond angle. In the magnetically unfavourable 
ferromagnetic case the energy of the system can therefore 
be lowered by distorting the bond angle away from 180°. 

It is clear from these results that qualitatively LaCr03 
behaves very similar to CaMnOs. However, on a quan- 
titative level, the polar instability is significantly weaker 
in LaCr0 3 than in CaMn0 3 (w 2 = -0.1 1 • 10 5 cm" 2 for 
+3% expansion in lattice constant, compared to uj 2 — 
-0.36 • 10 5 cm" 2 for +2% expansion in CaMn0 3 ). Even 
at rather large volume, the polar instability in LaCr03 
never dominates over the antiferrodistortive modest 



B. Born effective charges 

To further analyse the differences and similarities be- 
tween LaCr03 and CaMn03, we calculate the Born effec- 
tive charges (BECs) of the S-site cations in both systems. 
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TABLE I: Born effective charges Z* z of the TM cations in the 
three systems under consideration, calculated using the Berry 
phase approach, and compared to previously reported values 
and to the corresponding formal charges. 





BaTi0 3 


LaCr0 3 


CaMn0 3 


formal charge 


4 


3 


4 


U cB = eV 


7.31 


3.62 


7.61 


Uee = 4 eV 




3.58 


7.66 


previous work (U c s = eV) 


72 pa 


3.7rpl 


S.lfP 




6.89^ 




6.9P 



The BEC tensor Z* a a describes the change in polarisa- 
tion component P a resulting from a displacement of ion 
i along Cartesian direction /3! 2 ° l 21 l 



n dP a 

|e| dr it /3 



(1) 



Here, O is the unit cell volume and e is the electronic 
charge. Anomalously enhanced BECs (compared to the 
formal charge value) indicate that the corresponding dis- 
placements lead to strong changes in hybridisation with 
the surrounding ions, and are generally interpreted as 
"smoking gun" for the charge transf er mech anism to- 
wards ferroelectricity described in Sec. [I| 20 l 21 l 

We calculate BECs from the change in polarisation 
corresponding to small finite displacements of the TM 
cations along the cubic axes. The electric polarisation is 
calculated in two different ways: i) by using the Berry 
phase approach] 23 ! 24 ! and ii) by constructing maximally 
localised Wannier functions (MLWFs) and monitoring 
changes in t he cen tres of gravity of the individual Wan- 
nier 

orbitalsPEH The BECs obtained using the Berry- 
phase approach are listed in Table |TJ In addition to the 
BECs of the Mn and Cr cations in CaMn03 and LaCr03, 
respectively, the BEC for the Ti 4+ cation in the proto- 
typical ferroelectric BaTi03 is also calculated.^ 

It can be seen from Table[I]that the effect of U e g on the 
BECs is small and that all calculated values agree well 
with previously reported data. Furthermore, in agree- 
ment with Ref. El the BEC of Mn 4 + in CaMn0 3 is very 
strongly enhanced compared to the formal charge of +4. 
The anomalous part of the BEC amounts to about 90 % 
of the formal charge (or ~3.65 electrons), very similar to 
the enhancement observed for the Ti 4+ cation in BaTiOs. 
On the other hand the BEC of Cr 3+ in LaCr0 3 is only 
increased by a factor of 1.2 (or 0.6 electrons) compared 
to the formal charge of +3. Thus, the differences in 
the BECs between LaCr03 and CaMnOs seem consistent 
with the significantly weaker polar instability of LaCr03 
compared to CaMnOs. However, we also point out that 
there is no straightforward relationship between the mag- 
nitude of the BECs and the tendency of a certain mate- 
rial towards ferroelectricity (otherwise CaMnOs should 
be just as ferroelectric as BaTiOs). 



TABLE II: Decomposition of Born effective charges Z* z of 
the B site cation in contributions of individual Wannier cen- 
tres. Contributions of oxygen-centred MLWFs are averaged 
over both spin-projections, whereas TM ti g contributions cor- 
respond to the local majority spin-projection. All values are 
calculated for U e g = 0. 







BaTi0 3 


LaCr0 3 


CaMn0 3 


Mn/Cr 




— 


-1.081 


-1.019 


Mn/Cr 


d xz /d yz (f) 




-1.111 


-1.002 


0( 2 ) 


Px/Py 


0.853 


0.239 


0.516 


0(2) 


Pz 


0.330 


0.347 


1.067 


0(2) 


s 


0.267 


0.212 


0.331 


0(2/) 


Pi- 


-0.0812 


-0.059 


-0.088 


0(y) 


py 


-0.141 


-0.217 


-0.300 


0(2/) 


Pz 


- 0.092 


0.023 


0.165 


0(2/) 


S 


-0.072 


-0.129 


-0.124 


total valence 


3.065 


-2.760 


0.451 


semi-core 




-7.716 


-7.628 


-7.877 


ionic core 




12.000 


14.000 


15.000 


total BEC 




7.349 


3.612 


7.573 


formal 




4 


3 


4 



C. Wannier decomposition of BECs 

To further analyse the origin of the enhanced BECs in 
the three investigated systems, we decompose the electric 
polarisation calculated via MLWFs into contributions of 
the individual Wannier orbitals |u> n )P^ 



P -f'ioiiic core q ^ ^ (^rel^l^n 



(2) 



This translates into a corresponding decomposition of the 
BECs: 
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with: 



' Ar in 



(3) 



(4) 



Here, Ar„ = (w„(Ar ion )|r|u! n (Ar ion )) - (w n {0)\r\w n (0)) 
is the displacement of the Wannier centre n resulting 
from the ionic displacement Ar*i on . 

For each of the three systems under consideration, the 
occupied valence and semi-core states form energetically 
isolated groups of bands with a specific dominant atomic 
and orbital character: TM semi-core 3s and 3p, j4-site 
cation semi-core s and p, and oxygen 2s and 2p. In the 
case of LaCrOa and CaMnOs, the oxygen 2p bands are 
inter-mixed with the local majority spin TM t2 g states. 
We construct separate MLWFs for each isolated set of oc- 
cupied bands using the Wannier90 codeP^l The resulting 
MLWFs exhibit a clear atomic and angular momentum 




BaTiOo 



FIG. 2: Geometry for the calculation of BECs listed in Ta- 
ble [II] 0(y) and O(z) are the oxygen anions situated adjacent 
to the displaced TM cation along the y and z directions, re- 
spectively. The thick arrow indicates the displacement direc- 
tion (z direction). Picture generated using VESTAp^ 



(lm) character, albeit with a certain amount of admix- 
ture of other atomic orbitals on the surrounding ions due 
to hybridisation (see Fig. [3] for some examples). We use 
this atomic and orbital character to label each MLWF. 
The individual contributions to the TM cation BECs of 
all individual valence MLWFs as well as the sum over 
all semi-core contributions are listed in Table [II] 0(z) 
and O(y) denote the two symmetry inequivalent oxygen 
anions situated along and perpendicular to the displace- 
ment direction of the TM cation (see Fig. [2]). It can be 
seen that the values for the total BECs calculated via 
MLWFs which are listed in Table [n] are in good agree- 
ment with the values calculated via the Berry-phase ap- 
proach listed in Table [I] The individual contributions to 
the BECs for BaTi03 are also in good agreement with 
results of a similar decomposition presented in Ref. 1301 

Within a completely ionic picture, i.e. without inter- 
site hybridisation, the contributions to the BECs from 
the occupied TM ti g orbitals in LaCrOs/CaMnOa would 
be exactly equal to —1, whereas the contributions from 
all oxygen orbitals would be identically zero. In addi- 
tion, with the core/valence separation used in our pseu- 
dopotentials, the total contribution of all semi-core states 
would be exactly equal to —8 for all three systems. 
The t2 g contributions to the BECs in both LaCrOs and 
CaMnC>3 are indeed very close to the nominal value of 
— 1, and the semi-core contributions to Z* are also close 
to —8, deviating only by about 0.1-0.4 electrons from this 
value. In contrast, strong anomalous contributions are 
found for the oxygen s and p orbitals, with the largest 
contributions resulting from the p orbitals correspond- 
ing to the oxygen anions O(z) situated above and below 
the B-site cation along the displacement direction (z- 
direction, see Fig.[2|. 

The large positive contributions of these orbitals indi- 
cate that the centres of the corresponding MLWFs shift 
towards the TM cation that is moved closer to the oxy- 
gen. This represents a net electron flow from the oxygen 
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FIG. 3: (Color online) Maximally localised Wannier func- 
tions corresponding to 7r-oriented p x orbitals ((a), (b)) and a- 
orientedp z orbitals ((c), (d)) centred on O(z) for BaTiOs (left 
column), LaCrO.-; (middle column) and CaMnOs (right col- 
umn), (a) and (c) correspond to the cubic structures, whereas 
in (b) and (d) the TM cations were displaced along +z by 1% 
of the cubic lattice constants. Shown are cuts through the 
x-z plane. The positions of the oxygen anions/TM cations 
are indicated by the filled small/large (red/blue) circles. 



anion towards the TM cation. Fig. [3] depicts the changes 
in the O(z) p MLWFs for a displacement of the TM cation 
along +z by 1% of the cubic perovskite lattice constant. 
We point out that this displacement is much larger than 
the one used for the calculation of the BECs (which was 
only ~^0.2 %). A larger displacement was used in Fig.[3]to 
make the changes in the MLWFs more visible. It can be 
seen that the displacement of the TM cation towards the 
oxygen anion increases the amplitude of the correspond- 
ing atomic d character in the O(z) p MLWFs compared 
to the cubic case. This change in hybridisation shifts 
the centre of the Wannier function towards the closer 
TM cation, leading to the anomalous contribution to the 
BECs. 

In BaTiOa, where the Ti d states are empty and there- 
fore do not contribute to the polarisation, the anomalous 
contribution to the BECs listed in Table [TT] are strongest 
for the 7r-type (p x /Py) orbitals, i.e. for the p orbitals 
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that are oriented perpendicular to the Ti-0 bond and 
hybridise with the empty t^g states on the Ti cation. 
The contribution of the cr-type {p z ) orbital, which hy- 
bridises with the Ti e g states, is more than a factor of 
two smaller. 

If we compare this to the case of LaCrC>3, where the 
majority spin t 2g states are filled, the relative contribu- 
tions of the two types of p states change significantly. 
The contribution of the p z orbital is very similar to the 
case of BaTiC>3, but the contributions of the p x and p y 
orbit als are strongly reduced and are now smaller than 
the p z contribution. 

This effect is related to a strong reduction of the d-like 
"tails" in the O(z) p x and p y MLWFs in LaCrC>3 com- 
pared to BaTi03, which can be seen from Fig. [3^, and 
Fig. [3)3. In BaTiC>3, these tails represent the TM d char- 
acter contained in the nominal oxygen p bands, which, 
as discussed in the introductory paragraphs, are bond- 
ing states that result from hybridisation between atomic 
TM d and oxygen p orbitals. Once the corresponding 
antibonding states (or parts thereof) become occupied, 
the different orbital contributions can be separated into 
different MLWFs with appropriate orbital character. 

More specifically, we note that the MLWFs shown in 
Fig. [3] correspond to the global "spin- up" projection. For 
BaTiC>3 the two spin-projections are of course identi- 
cal, whereas for the magnetic systems the "spin-up" and 
"spin-down" MLWFs centred at O(z) are related to each 
other by space inversion. For the specific cases shown in 
Fig. [3] the local magnetic moment of the Cr/Mn cation 
located at +§£ relative to the central oxygen anion is 
parallel to the global spin-up direction, whereas the local 
magnetic moment of the Cr/Mn cation at —§2 is paral- 
lel to the global spin-down direction. It can be seen that 
the t2g-tails of the p x MLWFs at +|£ have completely 
vanished. This is due to the effect described above, i.e. 
the corresponding orbital character has been transferred 
to the t 2g {t) MLWFs. In addition, the t 2g -tail at -|i 
is also significantly reduced, because the global spin-up 
direction at this TM site correspond to the local minority 
spin character. The corresponding t 2g states are therefore 
higher in energy, which reduces the amount of p-t2 g hy- 
bridisation at this site. As a result, the p x /p y MLWFs in 
LaCr03 resemble more closely the corresponding atomic 
orbitals compared to the p x /p y MLWFs in BaTiC-3, and 
the contributions to the BECs become less anomalous. 

The same hierarchy between the p x /p y and p z contri- 
butions as in LaCrOs can be observed for CaMn03, but 
the overall magnitude of both contributions is strongly 
enhanced in the latter. The contribution of the O(z) p z 
MLWF in CaMn.03 is even three times larger than the 
corresponding contribution in BaTi03. 



D. Densities of states 

The origin of the strong enhancement of the individ- 
ual contributions to the BECs in CaMnC>3 can be ra- 
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FIG. 4: (Color online) Spin-resolved total and projected DOS 
fur cubic BaTi0 3 (top), LaCr0 3 (middle), and CaMn0 3 (bot- 
tom), calculated within GGA. The shaded curves represent 
the total DOS, the bright/dark (red/blue) lines represent the 
t2g/e. g projected DOS of the B-site cations. Different spin 
projections are shown with different sign. For better compar- 
ison, the DOS are aligned to the top of the oxygen p bands. 
Mid-gap (Fermi) levels are indicated by dashed vertical lines. 



tionalised by looking at the differences in the electronic 
structure of the three materials. The total and pro- 
jected densities of states (DOS) for BaTi03 (nonmag- 
netic), LaCr0 3 , and CaMn0 3 (both G-type AFM) are 
shown in Fig. [4] It can be seen that in BaTi03 and 
LaCr03 the TM d bands are situated well above the 
oxygen p bands. In particular the bottom of the TM 
e g states is about 3.5-4 eV above the top of the pre- 
dominantly oxygen p bands in both systems. On the 
other hand in CaMn0 3 the corresponding energy sepa- 
ration is much smaller, with the bottom of the Mn e g 
states only about 0.5 eV above the oxygen p bands. This 
small energy separation leads to very strong hybridis- 
ation between oxygen p and TM d levels in CaMn03. 
The much stronger hybridisation in CaMn03 compared 
to both LaCrOs and BaTi03 will thus enhance all result- 
ing effects, including the anomalous contributions to the 
BECs. The strong hybridisation in CaMn03 therefore 
compensates the reduced p x /p y contribution relative to 
p z , leading to the large total BEC (of the same magni- 
tude as in BaTiOs) and the rather strong polar instability 
compared to LaCrOs. 

As a final technical note, we point put that the weak 
effect of U e ff on the BECs in CaMn03 and LaCrOs is a re- 
sult of the specific projection on (non-orthogonal) atomic 
orbitals used for the LSDA+C/ implementation within 
Quantum ESPRESSO. For both LaCr0 3 and CaMn0 3 
the resulting orbital occupation matrix element corre- 
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sponding to the majority spin e g states is approximately 
equal to 0.5. Since the LSDA+[/ potential shift of the 
TM d states is given by AV£ = ^(§-rc£ l )P 1 the energy 
separation between the e g states and the oxygen p bands 
is thus not affected by the value of U e g for n| ~ 0.5. Ac- 
cording to the preceding discussion it is essentially this 
energy difference that determines the magnitude of the 
anomalous contribution to the TM BECs in LaCrC>3 and 
CaMn03, which are therefore rather independent of U e g. 

IV. SUMMARY AND CONCLUSIONS 

In summary, we have shown that LaCr03 exhibits a 
strongly volume-dependent polar instability very similar 
to CaMnOs. However the polar instability is significantly 
weaker in LaCrC>3 than in CaMn03. This is consistent 
with the much less anomalous BEC of the Cr 3+ cation 
in LaCrC>3 compared to the strongly enhanced BEC of 
Mn 4+ in CaMnOs. By decomposing the BECs in contri- 
butions of individual Wannicr functions, we could show 
that in both cases the largest anomalous contributions 
stem from er-oriented O p states of the oxygen anions 



adjacent to the TM cation along the displacement di- 
rection, whereas the corresponding tt contributions are 
reduced compared to the d° ferroelectric BaTi03. This 
reduction of the 7r-type contributions to the BECs rela- 
tive to the cx-type contribution is due to the occupation 
of the majority spin ti g states in LaCrC-3 and CaMnOs. 
The strong overall enhancement of the anomalous charges 
in CaMnC>3 can be explained by the near degeneracy of 
O p and TM d states in CaMn03, which leads to very 
strong hybridisation and thus enhances all related effects. 
The polar instabilities in both magnetic perovskites can 
therefore be understood in terms of charge transfer be- 
tween TM d and O p states, analogously to the standard 
d° perovskite ferroelectrics. 
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